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While all bipartite pure entangled states violate some Bell inequality, the relationship between 
entanglement and non-locality for mixed quantum states is not well understood. We introduce 
a simple and efficient algorithmic approach for the problem of constructing local hidden variable 
theories for quantum states. The method is based on constructing a so-called symmetric quasi- 
extension of the quantum state that gives rise to a local hidden variable model with a certain 
numbers of settings for the observers Alice and Bob. We use this method to analytically construct 
local hidden variable theories for any bound entangled state based on a real unextendible product 
basis (UPB) with two measurement settings for Alice and Bob. The problem can be approached by 
semi-definite programming and we present our numerical and analytical results for various classes 
of states. 



It was John Bell |l| who quantified how measurements 
on entangled quantum mechanical systems can invalidate 
local classical models of reality. His original inequal- 
ity has generated a field of research devoted to general 
Bell inequalities and experimentally observed violations 
of such inequalities. 

Perhaps surprisingly, the nature of the set of states 
that violate local realism is poorly understood, although 
it is known from the seminal work of Werner that not 
all entangled states violate a Bell inequality. Recent re- 
sults in quantum information theory have revealed the 
complex structure of the set of entangled states but have 
as yet shed little light on the relation between this struc- 
ture and violation of Bell inequalities. For example, it 
has been conjectured by Peres [[| that so called bound 
entangled states which satisfy the Peres-Horodecki "par- 
tial transposition" criterion |Q, do not violate any Bell 
inequalities. There are various results that support this 
conjecture both in the bipartite and multipartite case, 
see Refs. |B 0, §, 0, but none of the results is 
conclusive. What has been lacking in the literature so 
far is a systematic way of deciding whether a quantum 
state does or does not violate some Bell inequality. The 
difficulty is that the possible types of local measurements 
and the number of measurements that observers can per- 
form is in principle unbounded and the enumeration of 
Bell inequalities is computationally hard ATP . 

In this Letter we present the first systematic approach 
for constructing local hidden variable theories for quan- 
tum states, depending only on the number of local mea- 
surement settings for each observer. Our approach has 
yielded both numerically constructed local hidden vari- 
able theories for a variety of quantum states as well as an- 
alytical results for Werner states Q and a class of bound 
entangled states based on real UPBs [^2| . 

Before we can state our main result, we recapitulate 
the mathematics of local hidden variable (LHV) models 
and Bell inequalities for bipartite systems |32|. We re- 



fer the reader to Ref. (§ [Tl], [b| for some literature on 
the theoretical formulation of general Bell inequalities. 
Each of the observers, Alice and Bob, has a set of local 
measurements. Let i = 1, . . . , s a be the number of mea- 
surements for Alice and let each measurement have o a (i) 
outcomes. Let k = 1, . . . , Sb be the number of measure- 
ments for Bob and Ofc(fc) be the number of outcomes per 
measurement. The probability Pij,u denotes the prob- 
ability that Alice's ith measurement has outcome j and 
Bob's fcth measurement has outcome I. A local hidden 
variable model assumes the existence of a shared random 
variable between Alice and Bob that is used to locally 
generate a measurement outcome depending only on the 
choice of the local measurement (and not on the choice 
of the other, remote, measurement). The local hidden 
variable model generates the probability vector P with 
entries Pij t ki when it generates measurement outcomes in 
accordance with these probabilities. Mathematically 
one defines a convex set S(s a , Sb, o a , Ob) which is the set 
of probability vectors P that can be generated by LHV 
models. It is known that S is a polytope and that the 
extremal vectors B of S are vectors with 0, 1 entries [13|. 
For more information on polytopes, see for example [14|. 
These extremal vectors B correspond to the situation 
in which the outcomes of the measurements are deter- 
mined with certainty and can be labelled by 2 sets of 
indices m = (mi, . . . , m Sa ) where rrii = 1, . . . ,o a (i) and 
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In words, each extremal vector specifies a single outcome 
with probability one for each local measurement, inde- 
pendently of the measurement made by the other parties. 
An example of an extremal vector for s a = s& = 2 and 
o a = Ob = 3 is given in Tabic HI 



For a quantum mechanical system p in Tid 
probability Pij.ki is given by 
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Here {Efi > : £V E tj = I <i A } are the POVM elements 
for Alice's zth measurement and {Eg} are the POVM 
elements for Bob's kth measurement. There is a viola- 
tion of a Bell inequality if and only if Pijjd cannot be 
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TABLE I: An example of an extremal B- vector for 2 settings 
for Alice and Bob with 3 outcomes per setting. 

In this Letter we will prove the first necessary con- 
dition for a state to violate a Bell inequality depending 
only on the number of settings for Alice and Bob. We will 
explicitly construct a LHV model in a (s a = 2, Sb = 2) 
setting for any bound entangled state based on a real un- 
extendible product basis JO]. Then we will discuss nu- 
merical work that shows that many of the known bipar- 
tite bound entangled states cannot violate a Bell inequal- 
ity with two settings for Alice and Bob. Finally we will 
partially reproduce and extend some of Werner's original 
results by showing that it is possible to use our procedure 
to analytically construct LHV theories for Werner states. 
It is noteworthy to mention that our methods (Theorem 
1 and Theorem 2) straightforwardly generalize to mul- 
tipartite states, even though we have not explored this 
direction. 

We will connect violations of Bell inequalities to the 
existence of a symmetric (quasi-) extension of a quan- 
tum state [j33| . An extension of a quantum state p on, 
say, a system AB, is another quantum state defined on a 
system ABC such that when we trace over C we obtain 
the original quantum state p. We are interested in the 
situation where the system C = A^^'^fl®'' 1 " 1 ' and 
we will demand that the extension be invariant under all 
permutations of the s a copies of system A among each 
other and similarly invariant under any permutation of 
the B systems. It is clear that if the quantum state p is 
separable, i.e. p = £ i J^(|V't){V'*IX4 ® (I0t)(0«l)sj sucn 
an extension always exists: we just copy the individual 
product states onto the other spaces: 

= Y^tHimm) 9 " ® g&x&d®"- (2) 

i 

If the state p is a pure entangled state, then it is also 
clear that such a symmetric extension cannot exist. The 
symmetry requirement implies that the pure entangled 
state pab must equal pa' b where A' is another A-system, 
which is impossible. In popular terms we may say that 



pure entanglement is 'monogamous', B cannot be entan- 
gled with A and A' at the same time. In some sense what 
we show in this paper is that (1) a violation of a Bell in- 
equality indicates that the entanglement in the quantum 
state is 'monogamous' and (2) there are many mixed en- 
tangled states whose entanglement is not monogamous. 

Thus the existence of a symmetric extension can be 
viewed as a separability criterion (see Ref. [p| for a 
similar but stronger separability criterion where one de- 
mands that the symmetric extension has positive partial 
transposes). For considering Bell inequality violations 
we generalize our criterion slightly and ask whether a 
state has a symmetric quasi- extension H p which is not 
necessarily positive. In order to define this notion we 
need the definition of a multi-partite entanglement wit- 
ness, which is an entanglement witness which can de- 
tect any multi-partite entanglement in a state. It has 
the property that for all states ipi , . . . , ip Sa , 4>\ , . . . , 4> Sb , 
(ipi, ■ ■ ■ ,i>s a ><l>i) ■ ■ ■ As b \H p \ipi,. . . ,ip Sa) (f>i, . . . ,(f> Si ) > 
0. 

Definition [Symmetric Quasi- Extension]: Let ir : 
H® s -> H® s be a permutation of spaces H in H® s . We 
define 

Sym(p) = ^Y^P^- ( 3 ) 

We say that p on Ha <8> ~Hb has a (s a , s&)-symmetric 
quasi-extension when there exists a multi-partite en- 
tanglement witness H p on H® Sa ® H% Sb such that 
Tr ® (»„_!) 8(. k -i)ifp = p and H p = Sym A ®Sym B {H p ). 

A ' B 

The reason for considering such quasi-extensions is 
clear from the following theorems which are the main 
results of this Letter. 

Theorem 1: If p has a (s a , Sb)-symmetric quasi- 
extension then p does not violate a Bell inequality with 
(s a , Sb) settings. 

Before proving this theorem, it is important to note the 
generality of the result; it holds for all possible choices of 
measurements which includes POVM measurements with 
an unbounded number of measurement outcomes. We 
will show below that the quasi-extension of p effectively 
creates a LHV model for p when Alice and Bob have s a 
and Sb arbitrary measurements. 

Proof We prove our theorem by extracting an LHV 
model from the quasi-extension. The LHV model for p for 
{s a ,Sb) settings should reproduce the vector Pij t u{p) — 
Tr Efj <g> E^p for all possible choices of POVM measure- 
ments {E^j,E^}, as a convex combination of the ex- 
tremal B- vectors, i.e. 

Pa,M = £jw({^>£h}>p)£££> ( 4 ) 
1 1 1 1 1 

where p m ,n(-) > 0. If a symmetric quasi-extension exists 
for p then TtE$ <g> Eg p = Tr [E$j ® E% ® I) H p . Us- 
ing the definition of the B-vectors, the properties of the 
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POVMs (£V Eff = I dA B ), and the symmetry proper- 
ties of Hp it is not hard to verify that 



P ij , k i(p)=TrEf j ®E*p: 
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and similarly for 



E^. Since H p is a quasi-extension p m , n ({E^, E^}, p) = 
TrE4«)EfiJ p > 0, and we have obtained a LHV model. 
□ 

One way of looking at this result is the following [jl6| . 
If p has a symmetric extension p, then instead of mea- 
surement on p, Alice and Bob can do measurements on p. 
Due to the symmetry Alice can do the first measurement 
on the first Alice space and the second measurement on 
the second Alice space etc. But now these measurements 
are all commuting, and can be considered as one big mea- 
surement. But we know that when Alice and Bob each 
have only a single measurement a LHV model for their 
measurements exists and thus we have a LHV model for 
the measurements on p. With this picture in mind, it 
is not hard to understand the following strenghtening of 
our results (see also Ref. fl7|): 

Theorem 2: If p has a (1, s;,)-symmetric quasi- 
extension then p does not violate a Bell inequality with 
Sb settings for Bob and any number of settings for Alice. 

Remark The theorem also holds when Alice and Bob 
are interchanged. 

Proof. The intuition behind this theorem relies on the 
fact that there are no violations of Bell inequalities when 
one party has only one measurement setting, thus sug- 
gesting that it is unnecessary to extend to copies of Al- 
ice's space as well as Bob's. Here is the local hidden 
variable model that we construct from a quasi-extension 
Hp, on Ha 



>H% Sb . We set 



p m , a (({Ef j ,E*},p)) = 



n*f =1 (ivi# 



(Tr I a <8> E n H p ) s 



(6) 



Each p m ,n is nonnegative since H p is an entanglement 
witness. We can substitute this expression in Eq. (|J) and 
verify that we obtain the correct probabilities P%j t H (p) by 
using the definition of the B-vectors, the normalization 
of the POVMs, and the symmetry of H p as before. □ 

This method for constructing LHV theories may be 
implemented both numerically and analytically. Let us 
first show a simple analytic construction of a (2, 2)- 
symmetric extension for any bound entangled state based 
on a real unextendible product basis ]l2] |. Let Pbe = 
I — J2i l a ij bi)(a-i, h\ be the projector onto such a bound 
entangled state, where {|ai, h) — |a*,6*)} is the real 
unextendible product basis. Our (unnormalized) exten- 
sion will be \^>)a 2 A 1 ®\®)b x B 2 ~Y,i \ a h a i: b i, b i )A 2 A 1 B 1 B 2 

where 1^) = J2i ^ is evident that this extension 
has the desired symmetry property. It is not hard to 



verify that by tracing over the systems A2 and B2 we 
obtain P%e = Pbe- The existence of a symmetric (2, 2)- 
extension implies the existence of both (2,1) and (1,2) 
symmetric extensions for the state by tracing out copies 
of A or B, so any Bell inequality violation for this class of 
states must involve more than two measurement settings 
for both parties. 

We have implemented numerical tests for the condi- 
tions of these two theorems. Firstly we look for the ex- 
istence of a symmetric extension with H p > 0. If such 
an extension does not exist, there is still the possibility 
that some other kind of quasi-extension does exist. We 
have focussed on the existence of a decomposable entan- 
glement witness H p because in both these cases the nu- 
merical problem corresponds to a semi-definite program 
§J|. We label the partitions of Uf a <g> H% Sb into bi- 
partite systems by p and we denote partial transposition 
with respect to one of the two subsystems as T p . A de- 
composable entanglement witness may then be written 
as H p = P + J2 P Qp" where P > 0, Q p > for all p. (In 
fact, due to the symmetry it is only necessary to include 
partitions unrelated by permutations of copies of A or 
copies of B in the sum, as in Ref. (h|.) 

Semi-definite programs correspond to optimizations of 
linear functions on positive matrices subject to trace con- 
straints. They are convex optimizations and are partic- 
ularly tractable both analytically and numerically. We 
show how to numerically construct symmetric extensions, 
the decomposable quasi-extension case is very similar. 
The condition that the partial trace of H p is p is equiva- 
lent to requiring that Tr (X (g> T)H p = Tr Xp for all oper- 
ators X on Ha ®Hb- If we write X in terms of a basis 
{<7j} for the vector space of operators then by linearity 
it is enough to check that this trace constraint holds for 
each element of the basis. We will assume that the basis 
is orthogonal in the trace inner product Tr a iO~j = dij and 
that Co = Id A ® Id B I V oIaAb- The index i ranges from 
zero to (dAds) 2 — 1- Consider then this semi-definite 
program 

minimize Tr K, 

subject to Tr Sym^ ® Sym B {oi ® T)K = r.;, i > 0, 

K>0, 

where r» = Tvaip. If the optimum is less than or 
equal to one, then, by adding a multiple of the iden- 
tity to the optimal K, we obtain some K p that satisfies 
TrK p = 1 as well as the other constraints. If we de- 
fine Hp = Sym^ <g) Sym B (K p ) it is clear that H p is a 
(s a , Sh)-symmetric extension of p. Duality properties of 
semi-definite programs imply that an optimum greater 
than one precludes the existence of a (s Q , Sf,)-symmetric 
extension Q (see also the Appendix). 

We have implemented this semi-definite program us- 
ing SeDuMi |l9| for several examples of bound entangled 
states with c?^ = ds = 3. The results are summarized 
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in Table [0]. For (1,2) settings the extension code took 
l-3s to run on a 500 MHz Pentium 3 desktop with 500 
MB of RAM, while the quasi-extension code took 1.5-4s. 
The computation for both forms of the code as described 
above should scale roughly as d 2 A 5Sa+4 d 2 B 5sb+A |§. How- 
ever as s a and Sf, grow, the permutation symmetry of 
the extensions can be used to dramatically reduce the 
size of the problem by block diagonalizing the matri- 
ces Sym^ (g) Sym B (pi ® I) and removing repeated blocks 
p0| , pl[ . For a fixed dA,ds the computation will scale 
polynomially with s a and We implemented such a 
code in the case of (1, 3)-extensions which took 1.5-4s 
for a single state. 

The Choi-Horodecki (C-H) states considered in Ref. 
p2| depend on a parameter a and include separable 
(a e [2, 3]), bound entangled (a S (3, 4]) and nonpositive 
partial transpose states for a > 4. They turn out to have 
(2, l)-symmetric extensions well into the range for which 
the states are entangled. Over the range a € [4.34, 4.84] 
they have decomposable symmetric quasi-extensions but 
no symmetric extensions showing that the former prop- 
erty provides a strictly stronger sufficient condition for 
the existence of an LHV theory. 

On the other hand, we found that the two pa- 
rameter family of bound entangled states introduced 
by Horodecki and Lewenstein p|| do not have (2, 1)- 
symmetric extensions or quasi-extensions. Also, many of 
the states described by Brufi and Peres p4j do not appear 
to possess symmetric quasi-extensions. However, for sev- 
eral examples of these two kinds of states we searched 
numerically over measurement settings to look for vi- 
olations of extremal Bell inequalities for s a — Sb — 2 
and three outcome measurements, and also s a = st, = 3 
and two outcome measurements, without success. Note 
that states may have (s, l)-extensions and no (1, s)- 
extensions, we have performed both tests in all cases. Al- 
though this possibi lity does not affect the overall conclu- 
sion, the states of [g4| for example are sufficiently asym- 
metric with respect to swapping A and B that for s = 2 
we found examples having one kind of extension but not 
the other, as well as states having both kinds (these are 
examples with (2, 2)-extensions which implies this latter 
condition). 

We found that, although they have (2, l)-extensions, 
only a few of the general complex UPB states of |2[| 
have (3, l)-extensions and similarly the C-H states have 
(3, l)-extensions for a reduced range of values of a. We 
did not find examples of Brufi-Peres states having (3, 1)- 
extensions. 

Finally we considered Werner states Q defined in di- 
mensions d = dA = ds > 2 as pw — dTz^{I(d — $) + 
(<i$ — l)V) where V is the flip operator. Werner Q 
showed that for $ > — 1 + these states do not vio- 
late any Bell inequality with an arbitrary number s a ,Sb 
of von Neumann measurements (in Ref. p6| the author 
constructs LHV models for arbitrary POVM settings for 



C-H |2|:a € 
Complex UPB |B 
H-L 23 1 
Brufi-Peres 24 
Werner |2 
Werner d — 2, <E> > 



(2,1), (1,2) 

ext q-ext 



[2, 4.33] 
yes 
no 
few 
d > 3 
-1/2 



[2, 4.84] 

yes 

no 

few 
d > 3 
-1/2 



(3,1), (1,3) 
ext 



[2, 4.00] 
few 
no 
no 
d > 4 
-1/3 



TABLE II: Numerical results on the existence of symmetric 
extensions (ext) and decomposable quasi-extensions (q-ext) 
for (s a = 1, Sf, = 2), (s a = 2,s 6 = 1), (s a = 1,8b = 3), and 
(s a = 3, s b = 1). 



a more restricted range of <3?). We found that using sym- 
metry techniques similar to those in Ref. [^7j it is possi- 
ble to analytically solve the dual optimization problem 
to the semi-definite program described above, see the 
Appendix. The value of the optimum establishes that 
all Werner states have symmetric extensions so long as 
s a + Sb < d. Hence these states have LHV theories for all 
Bell experiments where the minimum number of settings 
s = min(s a ,Sfc) satisfies s + 1 < d. This result is more 
general than Werner's in the sense that, like in Ref. p6| , 
it holds for general POVM elements, ft is weaker in the 
sense that the number of settings is bounded by the di- 
mension of the space. Numerical and analytical results 
(see Table II and the Appendix) show that Werner states 
for d = 2 actually have symmetric (quasi-)extensions be- 
yond this analytically derived bound. 

Even though our method is the most powerful tool to 
date for constructing local hidden variable theories, we 
believe that it is unlikely that every LHV model can be 
constructed from a symmetric quasi-extension. Our work 
is only the starting point for a more thorough exploration 
of the existence of (quasi-)extensions for entangled quan- 
tum states. In particular, it is an intriguing and open 
question whether there exist entangled states that have 
(s a = l,Sb — > oo) quasi-extensions. In fact we heard 
that it has been proved that only separable states have 
(s a — 1, Sb — > oo) extensions ]23, |2{| |3(J. 
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Appendix 

In this Appendix we discuss the semidefinite program 
that constructs symmetric extensions in more detail and 
show that it is possible to construct symmetric extensions 
for many Werner states using our semi-definite program- 
ming approach. The ingredients of the argument are 
the semi-definite programming duality, the behavior of 
convex optimizations under symmetry and some simple 
group representation theory for the permutation group. 

The semi-definite program that attempts to construct 
a symmetric extension has a natural dual that will turn 
out to be simpler to solve. We first review, for com- 
pleteness, some standard descriptions and properties of 
semi-definite programs and their duals. 

Vandenberghe and Boyd | fl8| write the general semi- 
definite program as 

minimize c T x, 
subject to F Q + Y,i XiFi > 0, 

where c is a vector of length m, and (Fq, Fi) are n-by-n 
Hermitian matrices with i = 1, ...,m. The optimization 
variables form the vector x, also of length m. If there 
exists any vector x such that F(x) = F + ^ ji XiF i > 
the semi-definite program is said to be strictly feasible. 
The dual optimization, also a semi-definite program, may 
be written 

maximize —TtFqZ, 
subject to Z > 0, 

Tr F t Z = Cj, 

where the optimization variable is the n-by-n matrix 
Z. Again, the dual semi-definite program is said to be 
strictly feasible if there is an Z > satisfying the trace 
constraints. 

The most important relation between the primal and 
dual optimizations is that all allowed values of the primal 
objective function are greater than all allowed values of 
the dual objective function. Thus feasible points of the 
dual problem can be used to bound the optimum of any 
semi-definite program. This property results from simply 
evaluating the difference between the primal and dual 
objective functions for any feasible pair (x, Z) 

c T x + TiF Z = J^Tr XiFiZ + Tr F Q Z = Tr^(x)Z > 0. 

i 

The first equality results from the dual feasibility con- 
straints and the inequality holds since Tr AB is positive if 



A and B are positive and Hermitian. Typically, the dual 
optimizations are more closely related, this is captured in 
Theorem 3.1 of ]l8| . If both the primal and dual forms of 
a semi-definite program are strictly feasible, their optima 
are equal and achieved by some feasible pair (x opt , Z op \). 

Comparing the semi-definite programs above with the 
semi-definite program for constructing symmetric exten- 
sions we see that it corresponds to the dual form defined 
above with Fi = Sym^ ® Sym B (ai ® I), Cj = r% = Traip 
and F a = I dA (g) I dB (g) I. We have m = (cU^b) 2 - 1 
and n = d^dg. The sign of the objective function was 
changed for clarity. Recall that there is a symmetric ex- 
tension for p so long as the optimum is less than or equal 
to one. From now on we will write Sym' = Sym^igiSym^. 
To see that this optimization is strictly feasible consider 
any, not necessarily positive, matrix K' satisfying the 
trace constraints (K' exists since these constraints fix 
only a small number of the matrix elements of K 1 ). Then 
K = K' +nId A ®Id B ®f- is strictly positive so long as n > 
is greater than the magnitude of the largest negative 
eigenvalue of K'. K defined in this way still satisfies the 
trace constraints so the semi-dcfinitc program is strictly 
feasible. If we define the matrix X = Id A ® Id B + S x i a i> 
the dual form of our semi-definite program may be writ- 
ten 

maximize 1 — TiXp 

subject to Sym'(A ® I) > 0. 

The advantage of this semi-definite program from the 
point of view of analytical work is that the number of 
variables is very much smaller and we will see in the fol- 
lowing that it can be further reduced by symmetry meth- 
ods. Numerical implementations solve both optimiza- 
tions at once. To see that this program is also strictly 
feasible consider the point a;, = for all i, for which 
A = Id A ® Id B an d so Sym' (A g) I) > 0. From this we 
can conclude, by Theorem 3.1 of Q as described above, 
that if the maximum of this optimization is less than or 
equal to one (Tr X opt p > 0) there is a symmetric exten- 
sion for p. It is not necessary to construct the extension 
explicitly. 

We now use the symmetry of the Werner states to 
simplify this semi-definite program. The Werner states 
have the property that (U <8> U)pw(U (8 UY = pw for 
all unitary transformations U on TL ^7j. As a result 
the objective function is unchanged under the action 
X -> {U®U)X{U®Uy since Ti (U(g>U)X(U(g>U^ p w = 
TrX(U U^pwiU U) = TrXp w . Similarly if we 
choose an X such that the positivity constraints are sat- 
isfied then 

E/®( s »+ s i>)Sym'(A: <g> l)i/t®(«»+»») > o. 

On the other hand, JJ®( Sa+Sb ) commutes with all matrices 
that permute the tensor factors and so it commutes with 
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the symmetrization operation. As a result 

Sym'((f7 ® U)X(U <g> C/) f ® I) > 0. 

So the set of allowed matrices X is also invariant under 
the operation X — ► (Z7 (8>{/)X(f ®?7 )^. In fact since sums 
of positive matrices are positive convex combinations of 
matrices (U ®U)X(U '(g) for different {7 are also feasi- 
ble and also achieve the same value of the objective func- 
tion. As a result we may restrict our attention to matrices 
X such that f(U®U)X(U®U)UU = X. For any 

feasible X of lower symmetry there is some feasible sym- 
metrized matrix X such that Tr Xpw = Tr Xpyj. As dis- 
cussed in ]|^| all the Hermitian matrices of this form may 
be written X = xi!^ A <g> Id B + xV for some real xj and x. 
For the matrices X arising in our semi-definite program 
Xj = 1 so we are left with a single variable optimiza- 
tion. (Note that there are general methods along these 
lines for reducing the dimensionality of semi-definite pro- 
grams with symmetry |po[ The key property is the 
convexity of the optimization.) 

The objective function to be maximized for the Werner 
state is 1— Tr ATpw = — Since pw is separable for $ > 

we already know that an extension will exist for positive 
$ and can assume that $ is negative in the following. As 
a result we wish to find the maximum value of x for which 
the following matrix inequality holds 

If <■''+•>> +xSym / (y®I)>0. 

The eigenvalues of Sym' (X 1) are simply 1 + x\i where 
Aj are the eigenvalues of Sym(V(8>I). Supposing that the 
largest magnitude negative eigenvalue of Sym'(V <g> I) is 
— A m , the optimum is x op t = 1/A m . Since the optimiza- 
tion is now over a single variable, the semi-definite pro- 
gram essentially reduces to an eigenvalue problem. The 
optimum of the semi-definite program will be — a; op t$ and 
if this is less than or equal to one there is a symmetric 
extension for pw So pw has a symmetric extension over 
the range — l/x op t = — A m < $ < 1. Note that the sym- 
metrization is a completely positive map and maps posi- 
tive matrices X to positive matrices Sym'(X ® I). Since 
X = I + xV is positive for -1 < x < 1, Sym'(X <g> I) 
must also be positive over this range so x op t > 1- 

To proceed further we must evaluate the smallest 
eigenvalue of Sym'(V <£) I). The matrix V swaps the 
Hilbert spaces belonging to Alice and Bob and there- 
fore is an element of the representation of the group of 
permutations S Sa + Sb that is made up of permutations of 
the s a + Sb copies of TL. If "Kaj) is the matrix that swaps 
the i-th and j-th spaces, we have V — 7i"(i,2)- For clar- 
ity we will modify the order of the Hilbert spaces that 
we used above and imagine that odd-numbered spaces 

1 = 1, 3, s a + Sb — 1 are copies of Alice's system and 
that even-numbered spaces are copies of Bob's. 

As a preliminary, consider evaluating Sym(<7(g>J®^ s_1 ^) 

on TL® S . Clearly <r® J®^ -1 ^ is unaffected by any permu- 



tation that fixes the first object and so in some sense the 
only permutations that matter are the ones that swap 
the first Hilbert space with one of the others. More for- 
mally, we can divide the elements of the permutation 
group of s objects into cosets of the subgroup that fixes 
the first object. Any permutation can be realized by 
some permutation that fixes the first object followed by 
the transposition it{i.i) for some i that labels which coset 
the permutation is in. Thus if we write the elements of 
the subgroup that fixes the first element as ff we have 

Sym(iJ) = -^ttHtt^ 

7T 

Finally, a(^lf^ s ^ maps to itself under all n so Sym(<7(g> 

lf { ) = E 7r (i,j)( cr ® I)^(i,i)/s. Applying this obser- 
vation to the local symmetrization in the particular case 
where X = V = vr( 12 ) gives 

Sym'0/01) = — X)X)7r (ai _ li2J) . 

i=l j = l 

So the symmetrized matrix reduces to a linear combina- 
tion of transpositions of Hilbert spaces. 

The matrices tt form a reducible unitary representation 
of the symmetric group S Sa + Sb . There is a single unitary 
transformation that will block-diagonalize all the 7r with 
the blocks being, possibly repeated, irreducible represen- 
tations of the group (see for example Since this is a 
single unitary transformation applied to all tt it does not 
change the eigenvalues of linear combinations of different 
it. Such linear combinations will be positive if and only 
if all the blocks are positive and so we can restrict our 
attention to the eigenvalues of the blocks of Sym'^ ® I) 
in this basis. Now suppose that the alternating represen- 
tation occurs in this decomposition at least once. Since 
in this irreducible representation even permutations are 
represented by 1 and odd ones by —1, the value of this 
block of Sym'(V 01) is — 1. This implies that x op t < 1- 
Combined with the earlier lower bound on the optimum 
we have a; op t = 1 and pw has a symmetric extension for 
all allowed values — 1 < $ < 1. 

It is a standard result that the alternating representa- 
tion occurs in the decomposition into irreducibles of the 
representation of 5 Sa + Sb by permutations of tensor fac- 
tors of u® {Sa+Sb) so long as s a + s b < d flj. In fact 
it is easy to see that the alternating representation oc- 
curs if and only if there is a non-trivial completely an- 
tisymmetric subspace of j{® i - Sa+Sb \ This completes the 
argument. All Werner states have a (s a , s&)- symmetric 
extension and an LHV theory for all Bell experiments 
with (s a , Sb) settings if s a + Sb < d. 
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Note that the converse is not true, a symmetric exten- 
sion may exist even if this condition is not met. Con- 
side r s n — s* — 2 and d — 2. — In both those cases the 



[9. 



deco mposition includes the two dimensional irreducible 

representation of £4 [31] that can be generated by 



[10; 



7T(1,2) 



7T(3,4), 7T( 2 ,3) 




Any block corresponding to this irrep is then equal to 



r (l,2) 



r (M) 



T(2,3) 



F (3,4) 



)/4 = 



in 



[12 



[13 



Sin ne this has eigenvalues ( — 1/2, 1/2) we have x op { 



It is straightforward to cheek the blocks corresponding 

to the other irreducible representations to confirm that 
a^opt — 2. As a result we may say that Werner states in 
two dimensions have (2,2), (2, 1) and (1, 2) -extensions if 
— 1/2 < $ < 1. It is straightforward to confirm numer- 
ically that the same statement applies to decomposable 
(2, 2)-quasi-extensions. The irreducible representations 
of S n for large n are typically high dimensional. Note, 
however, that further block diagonalization is possible in 
principle since Sym'(V (£) I) is invariant under any per- 
mut ation of Alice's spaces or of Bob's. This means that 



[14 

[15 



[16 
[17 



[18. 

M 



eac h block of Svm' (V <Z)D is in the commutant algebra of 



the reducible representation of S Sa x S Sb that arises from 
restricting an irreducible representation of S Sa + Sb to a 
representation of S Sa x S Sb . This commutant algebra will 
typically have very much smaller block size. 
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